Absence of Hole Confinement in Transition Metal Oxides with Orbital Degeneracy 
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We investigate the spectral properties of a hole moving in a two-dimensional Hubbard model 
for strongly correlated t^ g electrons. Although superexchange interactions are Ising-like, a quasi- 
one-dimensional coherent hole motion arises due to effective three-site terms. This mechanism is 
fundamentally different from the hole motion via quantum fluctuations in the conventional spin 
model with SU(2) symmetry. The orbital model describes also propagation of a hole in some e g 
compounds, and we argue that orbital degeneracy alone does not lead to hole self-localization. 
[Published in Phys. Rev. Lett. 100, 066403 (2008).] 
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One of the fundamental problems in solid state physics 
consists in understanding the motion of an electron or 
hole coupled to the other degrees of freedom in a mate- 
rial. In many cases, the other degrees of freedom (spin, 
orbital or phonon excitations) can increase the mass of 
the carrier and possibly localize it. For example, the un- 
doped parent compound of high-T c cuprates is an antifcr- 
romagnetic (AF) Mott insulator due to electron-electron 
repulsion. A hole doped into it was at first thought to 
be localized, because its movement would disturb the AF 
background and thus cost energy Only two decades 
later it was found out that quantum spin fluctuations 
heal the background and lead to a coherent hole mo- 
tion 0, IH, see Figs. QJa-b). This shows, how important 
it is to critically assess any approximation used and to 
identify possible mechanisms of the coherent hole motion. 

Novel aspects of localization occur in orbital models. 
The superexchange (SE) is then no longer SU(2) symmet- 
ric, and the lower symmetry Q leads to anisotropy and 
often to frustrated interactions [H . The resulting variety 
of possible scenarios [f| make compounds with orbital 
degree of freedom at once very interesting for material 
science and challenging to theory. The most relevant sys- 
tems of strongly correlated orbitals arc (nearly) degener- 
ate e g or t2 g orbitals. While e g orbitals are of interest 
in colossal magnctorcsistance manganitcs, t 2g bands are 
relevant to, e.g., cubic titanates Q or vanadates [§]. 

Two e g orbitals describe the two-dimensional (2D) fer- 
romagneticTFM) planes of LaMn0 3 @, K 2 CuF 4 [lfl] or 
Cs 2 AgF4 The ground state of the orbital SE model 
has Ising-like alternating orbital (AO) order, where quan- 



tum fluctuations are largely suppressed (12|. Neverthe- 
less, a hole doped into the AO background finds a way 
to move in manganites, see Figs. QJc-d). Here, we in- 
vestigate hole motion in a FM plane with two active t 2g 
orbitals (similar to manganites, spin excitations could 
contribute for the AF coupling between the FM planes 
[HI]), where AO order is also Ising-like and interorbital 
hopping is excluded. This situation arises, when a crystal 
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FIG. 1: (color) Hole motion in: (a-b) the spin model with 
SU(2) symmetry, and (c-d) the e g orbital model. Defects in 
the AF background (a) caused by hole hopping are healed 
(b) by spin fluctuations, while for a hole in the AO state (c), 
interorbital hopping does not generate defects at all (d) . 



field splits t 2g orbitals in d 1 or d 2 systems. If one of the 
three t 2g orbitals is either empty (d 1 ) or fully occupied 
(cP), the remaining two can form the AO order, as e.g. in 
the planes of Sr 2 V04 [3] with possible weak FM order 
fill . In addition to t 2g compounds, this "t 2g model" also 
describes e g orbitals in the above mentioned fluorides, 
where a crystal field induces d z 2_ x 2 /d 2 2_ y 2-type AO or- 
der I 111 ], so both quantum fluctuations and interor- 
bital hopping are quenched and cannot generate coherent 
quasiparticle (QP) propagation Q, shown in Fig. [TJ 

In this Letter, we show that a hole doped in a state 
with alternating t 2g orbitals is not confined but finds a 
way to move coherently via three-site effective hopping 
terms arising from SE, i.e., even in a model with strictly 
nearest-neighbor hopping. For the present orbital model, 
long-range hopping is not expected to be important, be- 
cause it is straightforward to verify that: (i) the sec- 
ond neighbor hoppings flip the orbital flavor [16| , so they 
do not contribute to QP dispersion, while (ii) the third 
neighbor hoppings which conserve the orbital flavor are 
considerably smaller than the three-site terms for realis- 
tic parameters. These latter SE terms are often neglected 
171 ] , but here they play a central role and determine QP 
propagation. This finding contradicts naive expectations 
of absence of coherent hole motion for the Ising-like SE 
in the present t 2g orbital model. We investigate the spec- 
tral function by a combination of analytic and numerical 
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FIG. 2: (color) Artist's view of possible hole hopping pro- 
cesses in the orbital model (fJJ with AO order (t <C U). 
Dashed (red) circles indicate initial (final) positions of the 
hole. The dotted lines connecting them give the hole paths, 
while the shaded (yellow) rectangles mark electrons that 
changed orbital due to hole motion. When the hole moves 
by hopping t (a), it has to turn by 90° in each step and thus 
creates a zig-zag string of frustrated bonds. In contrast, in (b) 
the hole moves along a direction via effective three-site terms 
Q generated by intermediate states with two holes separated 
by a double occupancy (blue cross). Part (c) shows the hole 
moving in the other direction after frustrating three bonds. 



methods to arrive at unbiased conclusions. 

A FM plane with ti g AO order is described by inter- 
acting spinlcss fcrmions which undergo one-dimensional 
(ID) hopping with conserved orbital flavor: 
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Here a\ (b\) creates an electron with flavor a (b) that 
is allowed to move by hopping t along b (a) axis and 
cannot hop along a (b) axis in a cubic system jj| ; U gives 
the energy of a doubly occupied site. Apart from its 
applicability to t2 g and certain e g orbital systems, the ti g 
Hamiltonian is of high theoretical interest, because it 
presents — to our knowledge — the only possibility to 
derive Ising SE from a Hubbard- like model [l8|. In the 
regime of large onsite Coulomb repulsion U 3> t, where 
the undoped system is a Mott insulator, it reduces to 
an orbital t-J z model in a similar way as the Hubbard 
model does reduce to the spin t-J model [l9[. The model 
Hamiltonian, 7it-j~ — T~tt + T~Lj z , consists of: 



n * = -* E a & - 1 E ~ b t~ b j > 

{ij}\\b {ij}\\a 

n ->* = ^E(^-i^) ■ 



(2) 
(3) 



Here a\ = aj(l — n«,) and b\ = b\{l — n ia ) are the cre- 
ation operators in the restricted space without double 



occupancies, as in the spin t-J model [19j, and the sum 
in Eq. ([3]) includes each bond (ij) only once. The cor- 
responding density operators are n ia and h^, with the 
total onsite density hi = h ia + h^. T z = ^(h ia — h^,) 
stands for the zth component of the pseudospin opera- 
tor and J = 4t 2 /U is the SE energy. A hole inserted 
into the b orbital of the AO ground state can move by 
hopping t only along the bonds (ij) || a, i.e., only in one 
direction, see Fig. [2l This first step costs excitation en- 
ergy Ei = | J, and further steps build a string of orbital 



excitations with ever increasing energy [Fig. 2(a)| . In 
contrast to the spin Hubbard model with isotropic hop- 
ping, the severe restrictions on hole hopping remove all 
mechanisms of healing the defects in the AO state - 
not only the quantum fluctuations but even the Trug- 
man loop processes [ioj], which would lead to a coherent 
propagation, are here excluded. Consequently, the spec- 
tral function A(k, u>) for the t-J z model obtained from 
the self-consistent Born approximation (SCBA) is inde- 
pendent of momentum k and consists of a ladder spec- 
trum (not shown), with well separated peaks similar to 
Rcfs. Hi]. 

Surprisingly, the spectral function for the full 
orbital model ([1]) obtained within the variational cluster 
approach (VCA) [2l| exhibits a distinct coherent low- 
energy mode [22| [shown by a solid line in Fig. [3ja)] , 
with a ID dispersion which depends on the orbital. 

The VCA is a variational method based on exact diago- 
nalization (ED) combined with the self-energy functional 
approach [231 ] . and has its roots in perturbativc cluster 
approaches for the Hubbard model [2J|. The method is 
ideally suited to the present problem, because it com- 
bines unbiased solution of the full Hubbard-likc Hamil- 
tonian ([1]) on a small cluster (here 10 sites) with access to 
the thermodynamic limit. Spontaneous symmetry break- 
ing in ordered phases is incorporated by optimizing an 
appropriate 'fictitious' field, in our case a staggered or- 
bital field. Consequently, the resulting optimal state has 
almost perfect (Ising-like) AO order. To control our re- 
sults, we also performed ED on 4 x 4-sitc clusters and 
likewise found a ID dispersion with similar bandwidth. 

The above results suggest that the t-J z model with 
its dispersionless ladder spectrum cannot reproduce the 
spectral density of the full model (JTJ) . The puzzle is re- 
solved by noticing that the complete low-energy model 
in second order includes also three-site terms 
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where {imj} denotes three adjacent sites in a row || a 
(or a column || b) with m in the middle. This effective 
hopping term is also obtained perturbatively from Eq. 
([1]) by allowing one double occupancy next to the inserted 
hole, and is thus again of the order oc t 2 /U as the SE 
term @ . Figure |2(b)| illustrates how a b electron moves 
over an occupied a orbital and interchanges with the hole 
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FIG. 3: (color online) Spectral function A(k, o>) for b orbitals 
at U — 8t along the main directions of the 2D Brillouin zone 
obtained for: (a) the t2 g orbital Hubbard model {1} in the 
VCA with a 10-site cluster; (b) the effective model Ti. c a © 
in the SCBA. First moments of the VCA spectra calculated 
using the energy intervals [0.9t, 2t], [—t, 0.9t], and [—2.5t, —t], 
are shown by solid and dashed lines in (a). 

without affecting the AO order. This leads to ID free 
propagation on two sublattices: e£ = cos(2k x ) for 
the one with occupied b orbitals, and = iJcos(2fc. y ) 

for the other. If we include H^J and treat the effective 
strong-coupling model 

n oS = n t + Hj* + n£ } , (5) 

the low-energy QP state indeed becomes dispersive simi- 
lar to the VCA, both in ED for a 4 x 4 cluster (not shown) 
and in the SCBA, see Fig. [3jb). However, the free dis- 
persion oc J is strongly renormalized, see below. For a 
single hole, the QP weight is almost k-independent in the 
SCBA. In contrast, in the VCA for the orbital Hubbard 
model ([1]) the spectral weight of the QP state decreases 
strongly between k = (0,0) and k = (tt, 7t), similar to 
the decrease of the spectral weight in the spin Hubbard 
model 25]. This reflects the fundamental difference be- 
tween the Hilbert spaces of the full (Q]) and the effective 
strong-coupling model ([5]). 

One finds that not only the first, but also all subse- 
quent peaks have distinct ID dispersions in the SCBA 
[Fig-Etb)]. We attribute this to a ID propagation anal- 
ogous to the QP state, but started after an even number 
of t hopping steps have generated string defects in the 
AO state. The positions of all peaks evolve from the ex- 
citations of the t- J z ladder spectrum where the hole is 
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FIG. 4: (color online) QP properties for the t^g orbital (t2 g ) 
and for the conventional SU(2)-symmetric spin model [SU(2)] 
for increasing J/t: (a) bandwidth W/t, and (b) weight oqp 
averaged over the 2D Brillouin zone. Diamonds (orbital) and 
circles (spin) denote the VCA data for the Hubbard-like mod- 
els, lines give the SCBA results for the orbital t- J z model with 
three-site terms (© (solid), as well as the usual spin t-J model 
with (dotted) and without (dashed) three-site terms. 



confined, and the pscudogap between the QP and the sec- 
ond peak scales as (J/t) 2 / 3 , similar to the spin t- J z model 
Q. Such processes cannot be properly included within 
the VCA using a 10-site cluster, where too few momenta 
k are available. The dispersive features of Fig. EJb) are 
thus replaced by several maxima and the spectral weight 
is transferred between them when k changes. Remark- 
ably, the first moments calculated for these structures 
(within the relevant energy regimes) show again similar 
k-dcpendcncc [Fig. [3£a)] as that found in the SCBA. 

Figure [4] compares QP features of the present t 2g 
model, obtained within the VCA and the SCBA, with 
those for the spin Hubbard/i-J model. In the latter case, 
the bandwidth W at small J (large U) is, as expected, 
approximately linear in J [26[, both in the VCA (Hub- 
bard model), as well as in the SCBA (t-J model), with 
and without three-site terms. The t-J model with three- 
site terms gives the complete second-order perturbation 
result for the Hubbard model, and consequently its band- 
width agrees with the VCA data for somewhat larger J 
(smaller U) than the t-J model without three-site terms, 
see Fig.HJa). In the t 2g case, where there is no bandwidth 
~ J coming from the quantum fluctuations, the propa- 
gation via three-site terms leads to a bandwidth that is 
nearly quadratic in J. The additional power arises from 
the structure of the wavefunction [27| . which renormal- 
izes the bare three-site bandwidth ~ J/4 by an additional 
factor oc J related to the QP weight oqp . It is instructive 
to take a closer look at the contributions of the three-site 
terms to the bandwidth of the usual SU(2) symmetric t- 
J model: The difference between the SCBA results with 
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FIG. 5: (color online) Spectral function for the ID t^g orbital 
model obtained in the VCA for U — 8t. Excitations at u> ~ 
1.2t and u ~ — 1.7t correspond to a hole doped into a localized 
a orbital, while a hole doped into a mobile 6 orbital gives a 
free propagation with dispersion LUk = jJcosfc + c. 



[dotted line in Fig. |Ua)] and without (dashed) three-site 
terms is almost exactly twice the ti g bandwidth. This 
factor 2 is due to the reduction by 1/2 of the ti g SE ((3|) 
when compared to the spin model. 

In case of the QP weight oqp [Fig. [4|b)] , the three-site 
terms have almost no impact on the Hubbard model and 
the SCBA results for both models do not differ strongly. 
While the VCA gives a similar weight for t^g orbitals as 
the SCBA, the values for the spin Hubbard model are 
considerably larger. Since the SCBA has been shown 
to give the same results as quantum Monte Carlo data 
for the t-J model extrapolated to the thermodynamic 
limit [28| , we believe that the weight given by the VCA is 
affected by finite-size effects, which appear to be weaker 
in the more classical tig model. 

Finally, we discuss the ID chain along a axis [Eq. (TT]) 
without the first term], because it gives even clearer in- 
sight into the role played by the three-site terms. In 
one dimension, electrons of only one (b) flavor can hop, 
similar to e g electrons in a ID FM chain [2§|. Without 
three-site terms, one finds just three dispersionless peaks, 
one arising from a hole in the mobile orbital and two from 
the localized orbital [29|. With three-site terms, the hole 
doped into a mobile b orbital propagates, see Fig. Its 
bandwidth is exactly J, i.e., we observe the full three-site 
dispersion because hole motion is here not renormalized 
by string excitations. For a hole doped into a localized 
a orbital, t-J z and Hubbard-likc model give two identi- 
cal dispersionless features, which arise from hopping to 
the sites next to the initial site of the hole. One might 
expect the three-site terms to have some effect, because 
they allow the hole to move over the entire chain instead 
of being confined to just three sites. However, the energy 
gain due to this derealization is small, and the corre- 
sponding spectral feature is invisible in Fig. [5] For the 

are not seen 



same reason, processes displayed in Fig. 2(c) 



in the 2D spectra of Fig. [3] cither, where they might be 
expected to show up as higher-energy excitations with a 
dispersion complementary to that of the first peak. 

Summarizing, we have investigated hole motion in a 
background with ti g - and in some cases e g - AO or- 



der, and have found a coherent hole motion via three-site 
terms. This mechanism is fundamentally different from 
the ones established so far in spin systems (quantum fluc- 
tuations), or for e g electrons in manganitcs (intcrorbital 
hopping). It can also be distinguished from hole motion 
via direct longer-range hopping terms because it behaves 
differently under particle-hole transformation. Further- 
more, the present model provides a realistic case with 
purely classical Ising SE interactions. As one still finds 
coherent hole motion, we argue that the hole confinement 
and dispersionless ladder excitation spectrum of the t-J z 
model — while being attractive mathematical ideas - 
are never realized in transition metal oxides. 
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